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^H Abstract 

We give an interpretation of the temperature in de Sitter universe in terms of a dynamical Unruh effect 
(3JT)j associated with the Hubble sphere. As with the quantum noise perceived by a uniformly accelerated 

observer in static space-times, observers endowed with a proper motion can in principle detect the effect. 
In particular, we study a "Kodama observer" as a two-field Unruh detector for which we show the effect is 
approximately thermal. We also estimate the back-reaction of the emitted radiation and find trajectories 
associated with the Kodama vector field are stable. 
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1 Introduction 

In Rcf. [1], Kodama introduced a vector field which generalizes the notion of time-like Killing vectors to 
space-times with dynamical horizons 1 . Recently many authors have focused their attention on the Kodama 
vector field and found interesting results related to the measurement of the Hawking radiation emitted by 
different kinds of dynamical horizons [3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14] or time-like naked singularities [15]. 

Here we present a more "physical" process of measurement in which the detector has finite mass and 
extension [16], and interacts with a massless scalar field like in Ref. [17]. The background is a de Sitter 
universe and the horizon is provided by its exponential expansion. The Kodama observer does not follow a 
cosmic fluid geodesic, but remains at constant proper distance from the horizon, and is thus "accelerated" 
with respect to locally inertial observers. Consequently, it perceives the de Sitter vacuum as a thermal 
bath with a temperature associated with the surface gravity of the Hubble sphere, much like a uniformly 
accelerated detector in Minkowski space-time measures a temperature associated with its acceleration. As 
usual, such result is obtained from the transition probability for the detector to emit a (scalar) radiation 
quantum and increase its energy. Further, assuming finite detector's size and mass, we show that the Kodama 
trajectory remains stable during this process. 

We shall use the metric signature (H ) and implement units where c = G = kg = 1. 
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1 For a clarification of the differences between Kodama and Killing vectors, see Ref. [2]. For the ability of Kodama vectors 
to define a local time flow, see Ref. [3] 
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2 Kodama observer in de Sitter space 



We recall the de Sitter metric for an observer comoving with the cosmic fluid at constant r reads 

ds 2 =dt 2 -e 2Ht (dr 2 +r 2 dn) . (2.1) 

We are here interested in an observer at fixed areal position which moves along the trajectory 

re Ht =K, (2.2) 

with K constant. This observer stays at fixed distance from its cosmological horizon and its four- velocity is 
given by 

u" = ; , (2.3) 

VI - H 2 K 2 V ; 

where k 11 = (1, — H r, 0, 0) is the Kodama vector related to the surface gravity of the de Sitter horizon 

2irT H = ^" 1 | X=H _ 1 =H . (2.4) 

We also note the observer's four-acceleration is 

H 3 K 2 H 2 Ke- Ht 



l-H 2 K 2 ' 1 - H 2 K 2 



,0,0 , (2.5) 



with modulus 



M 4 K 2 

One can also introduce coordinates naturally associated with the Kodama detector, namely K from Eq. (2.2) 
and T, so that the metric reads 

ds 2 = (1 - H 2 K 2 ) dT 2 + 2 H K dT dK - dK 2 - K 2 dQ 2 , (2.7) 

which shows that the Kodama vector k — d/dT is a Killing vector. 

3 Particle detectors 

The concept of particle becomes rather involved once gravity and curved space-times are taken into ac- 
count [18]. In the following, we shall therefore employ the particle detector proposed by Unruh [19] and 
DeWitt [20], which can be described as a quantum mechanical particle with energy levels linearly coupled 
to a massless scalar field via a monopole moment operator. 

3.1 Point-like detectors 

Let us consider a massless scalar field <f) with Hamiltonian obeying the massless Klein-Gordon equation 
with dispersion relation Wk = |k|. The point-like detector is a quantum mechanical system with a set of 
energy eigenstates {| E% )} which moves along a prescribed classical trajectory t = t(r), x = x(r), where r 
is the detector's proper time. The detector is coupled to the scalar field <j) via the interaction Hamiltonian 

fli„ = Ai(#), (3.1) 

with A a small parameter and M(r) the detector's monopole moment operator. 

Suppose that at time To the detector and field are in the product state | 0,.Eo) = | 0)| E ), where | 0) 
is the scalar field vacuum and | Eq ) a detector state with energy E [21]. We want the probability that at 
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a later time n > To the detector is found in the state | E\ ) with energy E\ ^ E , regardless of the final 
state of the field <j). The answer can be readily obtained in the so-called interaction picture where states are 
evolved by the Schrodinger equation with the interaction Hamiltonian H m t- To first order in perturbation 
theory, the amplitude for the transition from | 0, E n ) at time r = To to | ip, E\ ) at time t = t\ is then 

(<p,E 1 \0,E Q ) = ( <p, E 1 |fexp (-i 1 dr ffi„t(r)^ \0,E Q ) 

~ -tX(E 1 \M(0)\E o ) r dre lT ^- E ^(^\^(r)\0) . (3.2) 

J T 

where T is the time-ordering operator. Squaring (3.2) and summing over all final field states, yields 

^2\(<fi, E i I 0,^o ) | 2 = A 2 |M(0)|£ )| 2 f 1 dr n dr' e-^- E °^-r ')( |0(r)^(r')| 0) . (3.3) 

^ J To J To 

The pre-factor A 2 |( E\ |M(0)| E )| 2 depends only on the peculiar details of the detector, whereas the response 
function 

i? T0 , Tl (A£) = r dr r dr' e- lAE ^ T '><0 |0(t)0(t')| 0) (3.4) 

Jtq Jtq 

is the same for all possible detectors. Here, we have set the energy gap AE = Ex — E n ^ for excitations 
or decay, respectively. ,/From now on, we will only consider the model-independent response function. 

Introducing new coordinates u := r, s := t — t' for r > r' and u := r', s := t' — r for t' > r, the response 
function can be re- written as 

Rto^AE) =2 J 1 du ° dsRc(e- lAEs {0 \&(u) $(u - s)\ 0)) , (3.5) 

having used (0 \4>(t')<P{t)\ 0) = (0 |</>(r)^(r')| 0)*, since is a self-adjoint operator. Eq. (3.5) can be 
differentiated with respect to t\ in order to obtain the transition rate 



Rt q . t {AE) =2^ ° dsRe(e- lAEs (0 \$(t)$(t-s)\ 0)) 



(3.6) 



where we set n = r. If the correlation function (0 \^>(t)4>(t — s) \ 0) is invariant under r-translations, 
Eq. (3.6) can be finally simplified to 

R T0T {AE)= lim / ° ds e- tAjEs W^ e (x(s),x(0)) . (3.7) 

^0+ J-(t-t ) 

where 

= <° \«AT))i«™ o) . |x(T) _ x(T , )l , v ; , ; 2 T) _ i(r , ) _,.p <3.8) 

is the positive frequency Wightman function. 

3.2 Unruh effect in de Sitter universe 

Let us apply the above construction to Kodama detectors in de Sitter space-time. Since Kodama observers 
are just the stationary de Sitter observers at constant distance from their cosmological horizon, we know 
that they will perceive a thermal bath at de Sitter temperature [22] 

T * = ■ (3-9) 
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We can confirm this expectation by noting that de Sitter space is conformally flat with conformal time 
rj = —H~ 1 e~ Ht . Then, provided 1 — K 2 H 2 > 0, the Wightman function (3.8) becomes 

1 H 2 1 

w ^ x ') = -^2 e H{t+v) (1 _ K 2 H 2 ){e~ m - e- m ' - ie) 2 ' (3J0) 

and, in the limit of r — > — oo, the detector transition rate (3.7) reads 

A/An . H 2 ,. f +oc J exp(-«A£;Vl -K 2 H 2 s) , oii , 

where we also used 

T = t^Jl-K 2 H 2 . (3.12) 
As long as if > (expanding universe) we can write the denominator as 

(e Hs ' 2 -e- Hs ' 2 -ie"} = 4 sinh 2 [if (s - ze)/2] . (3.13) 
This function has infinitely many double poles in the complex s-plane, namely for s = sj with 

Si = ^3, jeZ. (3.14) 

Since we are interested in the case of AE > 0, we can close the contour of integration in the lower half 
plane, summing over the residues of all the double poles in the lower complex s-plane, with the exception of 
the s = pole which has been slightly displaced above the integration path by the ie-prescription. The well 
known result turns out to be confirmed, namely 

, /A N 1 AEy/l-K 2 H 2 

R(AE) = - -, r , 3.15 

2tt exp ( 2^CTg )_i 

showing the presence of a cosmological Unruh effect with de Sitter temperature (3.9) red-shifted by the 
Tolman factor — K 2 H 2 , here appearing as a Doppler shift due to the proper motion of the detector. In 
fact, we remarked before that K is also the static coordinate of the detector relative to the static patch. 
Another intriguing formula can be written if we recall the acceleration (2.6), namely 



. (3.16) 



VI - H 2 K 2 2ir 

The observed temperature here appears due to a mixing of a pure Unruh effect (the acceleration term) plus 
a cosmological expansion term (the H term), and is the de Sitter version [23] of a formula discovered by 
Deser and Levine for detectors in anti-de Sitter space [24]. 

Alternatively, we can understand this effect as the transition from cosmological energy AE, conjugated 
to cosmic de Sitter time t, to the energy AE \f\ — K 2 H 2 as measured locally by Kodama's observers, with 
proper time r of Eq. (3.12). 



3.3 Extended detector and back-reaction 

We now repeat the previous analysis assuming the detector's size is not a priori negligible. In order to 
simplify the computation, we still assume spherical symmetry and the detector therefore only extends along 
one dimension. We will not display all the details but focus on the main differences. 

It is easy to see that a detector moving along a Kodama trajectory (2.2) in an expanding de Sitter 
universe, has the same dynamics of a particle which moves along the separatrix in the potential of an 
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inverted harmonic oscillator. ^From the equation of motion (2.2), we can therefore introduce the effective 
Lagrangian 



L iho = m(^ + H 2 r 2 



(3.17) 



where m is a parameter with mass dimension whose relation with physical quantities will be clarified later. 
Our detector is initially (at t = 0) represented by a Gaussian wave packet of size b peaked around r = K, 



exp 



V>(0,r) 



.HKmr (r — K) 2 



ft 



2b 2 



(3.18) 



which is then propagated to later times by the propagator obtained from the Lagrangian (3.17) (see, e.g., 
Ref. [16]), 



G(t,r;0,r') = 



i H m 



2 tt ft smh(H t) 



exp 



Hm [(r 2 - r' 2 ) cosh(tf t) - 2rV] 
2ft smh(Ht) 



(3.19) 



The complete expression of the detector's propagated wavefunction is rather cumbersome, however we notice 
that its square modulus reads 



exp 



2 b 2 H 2 m 2 {r~Ke- Ht ) 2 



b 4 H 2 m 2_ h 2 + + b 4 H 2 m 2) cosn (ff t) 



(3.20) 



which represents a Gaussian wave-packet peaked along the classical trajectory (2.2) employed in Sections 3.1- 
3.2 and spreading in time. The classical behavior is properly recovered in the limit ft — > followed by 
b — >• [16], in which the detector's wavefunction ip(r, t) reproduces the usual Dirac ^-function peaked on the 
classical trajectory. 

Now, in order to study the probability for the detector to absorb a scalar quantum and make a transition 
between two different trajectories (parameterized by different rrii and if,), one needs to compute the tran- 
sition amplitude for finite b and ft (otherwise the result would automatically vanish, the response function 
involving the product of two Dirac <5's peaked on different trajectories). The detector also interacts with the 
quantized scalar field ip = <p(t, r) according to 



Ant = ^ Q (^2 V'l + "02 Ipl) f , 



(3.21) 



where Q is a coupling constant and ipi = ipi(t, r) two possible states of the detector corresponding to different 
trajectories Vj = Kie~ Ht and mass parameters rm 2 . We assume the difference between the two states is 
small, 



K X =K -\5K 
K 2 = K + \8K 



mi = m — \ 8m 
m 2 = m+ i 5m 



(3.22) 



and expand to leading order in 5K and 8m and, subsequently, for short times (Ht ~ Ht' <C 1), keeping ft 
and b finite. In particular, one obtains 



H 2 K 2 

-i—^—Sm (t-t') + Q{b) 



(3.23) 



2 A fundamental difference with respect to the Unruh effect analyzed in Ref. [17] is that the acceleration parameter H is not 
varied here, since it is a property of the background space-time. A change SK implies a change in the detector's acceleration 
according to Eq. (2.6). 
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in which we have evaluated the phase (in the saddle-point approximation) with r along the average trajectory 
between r\ and [17]. Like the Unruh detector in Rcf. [17], the above exponential does not contain a real 
(quadratic in 5K) part, contrary to the case of a geodetic observer [see Eq. (A. 5) in Appendix A)], which 
implies that the transition amplitude will not vanish in the point-like limit 6 — >• 0. Upon comparing with the 
result obtained for the point-like case, we immediately recognize that 



H 2 K 2 m, = E y/l - H 2 K 2 , (3.24) 
where E is the detector's proper energy and 



iplMt) MM*') - ex P 



- l -8E ^l-H 2 K 2 (t -t')+ l - 2 - R2 ; ESKtt- t') + 0(6) 



(3.25) 



We can now take the limit b — > 0, as part of the point-like limit in which one would not expect the second 
term in the above exponential. In Ref. [17], we required the analogue of the second term above vanished 
and obtained the equation of motion for a uniformly accelerated detector in Minkowski space-time, namely 
ma = f and constant. Following the same line of reasoning, we now obtain the equation of motion 

SK = 0. (3.26) 

This can be interpreted as meaning the Kodama trajectory is stable against thermal emission of scalar 
quanta in the de Sitter background. 

The transition probability per unit de Sitter time can finally be computed by taking the classical limit, 
in which one recovers the same result in Eq. (3.15) with AE — SE. 



4 Conclusions 

We considered an expanding de Sitter universe that is inflating exponentially and studied the dynamics of an 
observer as an object which is collapsing towards the Hubble sphere. This is the "Kodama observer" placed 
at K — re Ht constant, and its radial velocity (2.3) is in fact negative which means that the observer is 
moving towards decreasing radii. In particular, the observer's motion is described by an inverted harmonic 
potential and involves the negative radial acceleration in Eq. (2.5). We can picture our observer/detector in 
spherical coordinates as a "shell" with a "density" profile \ip\ 2 peaked on the average radius r [see Eq. (3.20)] 
which sees the universe becoming less dense (in time) but in a homogenous way (in space). Of course, the 
detector's proper time r in Eq. (3.12) differs from the cosmological time t for an observer comoving with the 
cosmic fluid. 

The detector interacts with a scalar field so that it can absorb or emit quanta and change its proper mass 
(energy) . Indeed, the Unruh effect is given by the simultaneous emission of a scalar quantum and detector's 
excitation. We therefore considered the transition between two states of the detector corresponding to two 
different trajectories associated with different energies and values of K and obtained the transition rate 
(3.15) in the point-like limit. Our result shows the expected thermal behavior (3.9) associated with the 
trajectory's stability (3.26). 

In summary, we studied the response of a Kodama detector moving along trajectories at constant distance 
from the Hubble sphere of the de Sitter universe and found that it perceives a thermal noise associated with 
the emission of scalar quanta. We also estimated the back-reaction of the emitted radiation and showed 
trajectories associated with the Kodama vector fields are stable. This represents a novel semiclassical 
property of the (classically defined) Kodama vector. 



A Extended detectors in flat space-time 

In order to clarify the analysis of the extended Kodama detector in de Sitter space, let us review the case 
of a detector moving along a geodesic in flat space-time. Without loss of generality, we consider a detector 
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at fixed x = in a Minkowskian coordinate system {t, x, y, z} which, for simplicity, extends only along 
the axis x. This case can also be interpreted as a detector comoving with the cosmic fluid in a general 
Friedmann-Robertson- Walker metric by replacing x with the radial comoving coordinate r. 
The initial wavefunction for such a detector is now 



cxp 



2F 



m 



2nih(t-t') 



and is propagated in time by the free particle propagator 
G(t,x;0,x') 

which yields 

ip(t,x) 



exp 



2 ) 



2h(t-t') 



i ^pii (ht — ib 2 m) 



cxp 



2(b 2 m + iht) 



(A.1) 



(A.2) 



(A.3) 



As before, we now consider the transition between two detectors with slightly different positions and mass, 



Xi 



X2 



5x 



h Sx 



mi 



m 2 



m 



\ Sm 



m + | Sm 



(A.4) 



and expand to leading order in Sx and Sm, keeping ft and b finite. In particular, by expanding in ft for b 
fixed, one obtains 



S_J_ 

4b 2 



ft t Sx 2 Sm 



rn. 



o(h 2 ) 



(A.5) 



in which we have evaluated the phase for x along the average trajectory between x\ and X2 (that is, x = 0). 
Note that the first term in the exponential generates a Dirac-<5(<5a;) which therefore forces the transition 
to occur at Sx — as one would expect for a geodetic motion with no emission or absorption. In fact, 
the condition Sm = is ensured by the Wightman's function evaluated along geodesies. Further, no terms 
proportional to 1/ft appear at all and, in particular, none analogous to the second one in Eq. (3.25). 
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